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ABSTRACT 

This paper investigates conditions under which balanced incomplete block designs enjoy weighted optimality 
with E-criterion establishing weighted intervals for E-optimal design. More so, a neighbourhood of weights for grouped 
generalised divisible designs (GGDDs) maintaining E w - optimal in D (v, b, k) was also investigated. The E-criterion was 
shown to be closely related to efficiency balance. Bounding arguments that are important tools in tackling E-optimality 
problems was employed; the standard bounds were generalized for seeking E-weighted optimal (E w -optimal) designs. 
The optimal bound established the best conceivable values of the criterion and thus the designs with these best values are 
optimal. 
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INTRODUCTION 

In the design of experiments, optimal designs are a class of experimental designs that are optimal with respect to 
some statistical criterion. According to Ipinyomi, (2012), under the design of experiments for estimating statistical models, 
optimal designs allow parameters to be estimated without bias and with minimum variance. A non - optimal design 
requires a greater number of experimental runs to estimate the parameter with the same precision as an optimal design. 
In practical terms, optimal experiments can reduce the costs of experimentation. The optimality of a design depends on the 
statistical model and is assessed with respect to a statistical criterion, which is related to the variance matrix of the 
estimator. Specifying an appropriate model and specifying a suitable criterion function both require understanding of 
statistical theory and practical knowledge with designing experiments. 

Gupta et al (1999) and Gupta et al (2002) used the term "weighted optimality" when comparing a group of test 
treatments with a group of control treatments. Two different sets of control, treatment - control and treatment - treatment 
were considered to be estimated with unequal precision. Design of experiments for which some of the treatments are 
controls has form a special perspective, been extensively investigated in recent years. Notable among the many papers 
seeking optimal designs for test treatment versus control experiments are Jacroux (1989), Jacroux and Majundar (1989), 
Majundar (1992, 1996), Majundar and Notz (1983) and Stutken (1991). Optimality work for T v C experiment is a limiting 
case, as the weight on the control treatment goes to 1 . 

Kiefer (1975) introduced convex optimality function § on the information matrices and proved that balance 
incomplete block designs (BIBDs) are universally optimal, i.e. minimise (j) (Cd) for every non - increasing, convex 
permutation - invariant §. Following closely on the heels of Kiefer's work, John and Mitchell (1977), Cheng (1978, 1980) 
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and Jacroux (1980) used computer search and theoretical arguments to build optimal designs for the criteria defined above, 
all of which fall into Kiefer's frame work. More recently, Majundar and Notz (1983), Majundar (1986), Jacroux and 
Majundar (1989), Bagchi and Shah (1989), Bagchi and Berkum (1991), Bagchi (1996), Bagchi and Bagchi (2001), 
Reck and Morgan (2005) and Morgan (2000, 2003, 2007) have been working on design optimality for various classes of 
designs with blocking. 

Let Y U j be the observation on experimental unit u in block j. The commonly employed statistical model for any 
block design d, which in many cases is justifiable by randomization alone (Hinkelmann and Kempthorne, 2008) is 

Yuj = /U + T d(u . j) +^ J + e UJ , (1) 

Where 

ji = mean response over all treatment and blocks, 

d( U ,j) = the treatment assigned to unit u in block j by design d, 

Td(uj) = the effect of the treatment assigned to unit u in block j by design d, 

/?j = the effect of block j, 

and the e a] 's are uncorrected, mean zero random variables with common variance er E 2 . This model is employed in 
most of the papers cited above. It is the basis for the information matrix mentioned earlier on. We usually assumed with no 
loss of generality that the unit variability cr E 2 is cr E 2 = 1 . The symbol n is used for the total number of experimental units, 
n = bk 

Weighted Optimality of Block Designs 

Kiefer's design optimality is based on functions of the information matrix that are invariant to treatment 
permutation, that is § (PC d P T ) = § (C d ) for any permutation matrix P. This implies equal interest in all treatments. 
However, in practice there are many cases where not all treatments are equally important. For instance, we often encounter 
experimental situations where some test treatments are to be compared to a standard treatment (or control treatment). 
Sometimes the control is included specifically to verify the expectation of large treatment effects relative to control, after 
which the important comparisons among test treatments are performed. This indicates asymmetry of interest on test 
treatments and the control treatment with (in this case) greater interest in test treatments than control. Asymmetry of 
treatment interest implies that optimality based on the information matrix should not be invariant to all permutations. 
According to Wang (2009), with the premise of asymmetric interest, the approach here is to group treatments into several 
subsets which are assigned distinct weights; larger weight reflects greater interest placed on estimating comparisons 
involving the corresponding treatments. In situations like that described above, this leads to a 2 - weight design problem, 
that is, the weights take only two values, with one small weight and v-1 larger, equal weights. 

Definition 1 

Let positive weights Wi, w 2 , . . . w v be measures of interest on v treatments without loss of generality 
Yu=i Wj = 1. Let W be a v x v diagonal matrix with w ; in the i th diagonal position, that is, W = Diag (w ; ). Also, the square 
root matrix for W is denoted as W' /2 . Then the weighted information matrix C dw for design d is defined as 
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c dw = w I/2 c d w I/2 = (( Cdii y Vw, w,o) (2) 

The use of C dw will be justified in the following facts. The key is to see how applying weights to C d induces 
weights on variances of treatment contrasts. Consider the spectral decomposition of the C dw -matrix: 

c dw = w ,/2 c d w I/2 = YX=l Uf i- (3) 

Where 9 < 9[ < . . .< 9 v _i are the eigenvalues of C dw and the fi are an orthonormal set of eigenvectors. 
For connected designs, both C d and C dw are of rank v - 1, and 9 = 0. Let w = (w b w 2i . . . , w v )' be the vector of weights. 
The eigenvector corresponding to 9 is/ = w' /2 = (V w b V w 2 , . . . , V w v ,)' 

Basics for the Weighted E- Criterion 
Definition 2 

The weighted E - value (written as E w ) for a design d is the largest canonical weighted variance for design d. 

That is, 

E w =± (4) 

A design d is weighted E'-optimal (or ^-optimal) in a design class D if it minimises the largest canonical 
weighted variance, that is, if 

E- dw = min£ dw (5) 
dED 

Result 2: The weighted fi-value is the largest weighted variance over all treatment contrasts. 
Proof 2: The largest weighted variance over all contrasts is 

max ( ) = max ( — ) = max ( — ) 

cWc cW c y'y 

c'l=0 c'l=0 y'w il2 =0 (6) 

w 1 ' 2 is an eigenvector of Cj" w corresponding to eigenvalue 0, so this is the largest eigenvalue of Cj" w , that is l/9i. 

Result says that an E w -optimal design factors importance of contrasts into design selection in minimising impact 
of the worst case. It can be seen that for any design, increasing the weight placed on a treatment increases weighted 
variances Var lv (c'x) = [c'V^'c] -1 Var (c't) of contrasts in which it is involved. Minimising summary functions of weighted 
variances (that is, minimising functions of 1/9,), pushes variances of treatments with higher weight to be smaller, this being 
at the expense of variances of treatments with smaller weights. 

STATEMENT OF THE PROBLEM 

The general formulation on the studies of optimality for treatment comparison is based on the idea that optimality 
functions of the treatment information matrix are invariant to treatment permutation which implies equal interest in all 
treatments. In practice, however, there are many experiments where not all treatments are equally important. 
When selecting a design for such an experiment, it would be better to weigh the information gathered on different 
treatments according to their relative importance and/or interest. It is on this premise that this research work is based on 
weighted optimality criteria. 
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AIM AND OBJECTIVES 

The broad objective of this research is to explore weighted optimality of incomplete block designs. The specific 
objectives are: 

• To establish weighted intervals for E -weighted optimal (E w - optimal) design covering possible Treatment with 
Control (T W C) situation having smaller weight on contrast. 

• A neighbourhood of weights for Grouped Generalised Divisible Designs (GGDD) maintaining E w - optimality in 
D (v, b, k) is also investigated. 

MATERIALS AND METHODS 

The conventional criteria for evaluating design optimality are functions of the eigenvalues of the information 
matrix C d , in the same way, many of the weighted criteria which is used to evaluate design optimality are function of the 
weighted information matrix C dw . There are many statistical packages for analysing the information matrix of a design 
optimality such as MINITAB, SPSS and R - Statistics. In this research, I made use of R - Statistical package. 

ANALYSIS OF RESULTS 

This chapter deals with the analysis and interpretation of data. I considered different designs of different sizes for 
E w -optimality criteria and a neighbourhood of weight for grouped generalised divisible designs (GGDD). 

Analysis and Results of Optimal Block Designs 

Example 1: For Vi = 6, v 2 = 3, Wi = V 12 and w 2 = the following design is E w - optimality over 
D(v, b, k) = (9, 11,6) 
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It can be observed that the above design is built up by adding the treatment 7, 8 and 9 to every block in 
BIBD (6, 10, 3), then appending symbols (1, . . . , 6). 

Example 2: The following design in D (7, 17, 4) is E w - optimal for Vi = {1, 2, 3}, v 2 = {4, 5, 6, 7}, Wi = Vi and 

w 2 = V 16 

llll 1 111111111114 
222 2 2 222222222225 
3 3 3 3 3 33333333333 6 
4444 : S 5 5 6 6 6 6 7 7 7 7 7 



Grouped Generalised Divisible Design 
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Example 3: The following two designs are built up from BIBD (7, 7, 3). The first design adds b = 1 copy of first 
block in BIBD (7, 7, 3) to obtain a GGDD (2). The second one adds b = 2 copies of the first block in BIBD (7, 7, 3). 

For w t = 2, w 2 = 1, Vi = 3 and v 2 = 4, this design is E w - optimal in D (v, b, k) = (7, 8, 3). 



1112 2 3 3 1 

24643452 
3 5 767763 



For Wi = ¥i3, w 2 = Vi3, Vi = 3 and v 2 = 4, this design is E w - optimal in D (v, b, k) 



1 1 

2 4 
a 5 



Example 4: Consider the following BIBD (5, 10, 3) 
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Y, = Y5, observed that 2 / 3 > k / v > V 2 , so a = 1 (by theorem). It can be checked that X (v — k - 1) < 2 (k - 1), so d can 
be any positive number. The designs constructed by appending b copies of one block in the above BIBD are E w - optimal 
in D (v, b, k) = (5, 10 +b, 3) for Vl = 3, v 2 = 2 and w 7w 2 = (X +b)fk = b / 3 + 1. 

Table 1: Parameters of Weight Balanced, Binary Block Designs 
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Table l:Contd.. 
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CONCLUSIONS 

Weighted intervals for Ew - optimality was established and it was shown to be closely related to efficiency 
balance. More so, bounding arguments establishes the best conceivable value of E-criterion in the first main result for 
Ew-optimal ruling out classes of inferior designs and since each of these designs achieved the best conceivable value 
respectively and the E w - optimal designs for the 2 - weight problem with group sizes > 2 i.e. GGDD(2), the necessary and 
sufficient condition for d to have weighted information matrix Cd W = G (I - w 1/2 w 1/2 ') for some G and some w t and w 2 are 
met. Hence the designs considered in this paper are good designs. 

REFERENCES 

1. R. A. Ipinyomi, Lecture notes on design of experiments. 2012. 

2. V. K. Gupta, D. V. V. Ramana, & R. Parsad, Weighted A-optimal block designs for making treatment-control and 
treatment-treatment comparisons. Journal of statistical planning and inference!! , 1999 301 - 319. 

3. V. K. Gupta, D. V. V. Ramana, & R. Parsad, Weighted A-optimal block designs for comparing test treatment with 
controls with unequal precision. Journal of statistical planning and inference 106, 2002, 159 - 175. 

4. M. Jacroux, the A-optimality of block designs for comparing test treatments with a control. Journal of the 
American statistical association. 84, 1989, 310 - 317. 

5. M. Jacroux, & D. Majundar, Optimal block designs for comparing test treatments with a control when k > v. 
Journal of statistical planning and inference. 23, 1989, 381 - 396. 

6. D. Majundar, Optimal designs for comparing test treatment with a control utilizing prior information. 
The annals of statistics. 20, 1992, 216 - 237. 

7. D. Majundar, Designs and analysis of experiments. Optimal and efficient treatment-control designs. 
The annals of statistics. 13, 1996a, 1007 - 1053. 

8. D. Majundar, on admissibility and optimality of treatment-control designs. The annals of statistics. 24, 1996b, 
2097 - 2107. 

9. D. Majundar, & W. Notz, Optimal incomplete block designs for a comparing treatments with a control. 
Annals of statistics. 11, 1983, 258 - 266. 



Index Copernicus Value: 3.0 - Articles can be sent to editor@impactjournals.us 



Weighted E - Optimality Criteria for Comparison of Balance Incomplete Block Design 



95 



10. J. Stufken, on group divisible treatment designs for comparing test treatments with a standard treatment in blocks 
of size 3. Journal of statistical planning and inference. 28, 1991a, 205 - 221. 

11. J. Stufken, Bayesian A-optimal and efficient block designs for comparing test treatments with a standard 
treatment. Communications in statistics. -.Theory and methods. 20, 1991b, 3849 - 3862. 

12. J. Kiefer, Construction and optimality of generalised youden designs. In Srivistava, J. N. A survey of statistical 
design and linear models. (North - Holland, Amsterdam, 1975), 333 - 353. 

13. J. A. John & T. J. Mitchell, Optimal incomplete block designs. Journal of the royal statistical society. 
Series B. Methodological. 39, 1977, 39-43. 

14. C. S. Cheng, Optimal designs for the elimination of multi - way heterogeneity. The annals of statistics. 6, 1978, 
1262- 1272. 

15. C. S. Cheng, On the E-optimality of some Block designs. Journal of the royal statistical society. 
Series B. methodological 42, 1980, 199-204. 

16. M. Jacroux, on the determination and construction of E-optimal block designs with unequal numbers of replicates. 
Biometrika 67, 1980a, 661 - 667. 

17. M. Jacroux, On the E-optimality of regular graph designs. Journal of the royal statistical society. 
Series B. Methodological.42, 1980b, 205 - 209. 

18. D. Majundar, Optimal designs for comparisons between two sets of treatments. Journal of statistical planning and 
inference. 14, 1986, 359 - 372. 

19. S. Bagchi, & K. R. Shah, on the optimality of a class of row-column designs. Journal of statistical planning and 
inference. 23, 1989, 397-402. 

20. S. Bagchi & Van Berkum, on the optimality of a new class of adjusted orthorgonal designs. Journal of statistical 
planning and inference. 28, 1991, 61 - 65. 

21. S. Bagchi, An infinite series of adjusted orthogonal designs with replication two. Statistical sinica. 6, 1996, 
975 -987. 

22. B. Bagchi & S. Bagchi, Optimality of partial geometric designs. The annals of statistics. 29, 2001, 577 - 594. 

23. B. Reck & J. P. Morgan, Optimal design in irregular BIBD settings. Journal of statistical planning and inference. 
129, 2005, 59 - 84. 

24. J. P. Morgan, Optimal designs with many blocking factors. Annals of statistics. 28, 2000, 367 - 387. 

25. J. P. Morgan, Optimal row - column design for two treatments. Journal of statistical planning and inference. 115, 
2003, 603 - 622. 

26. J. P. Morgan, Optimal incomplete block designs. Journal of the American statistical association. 102, 2007, 
655-663. 



Impact Factor(JCC): 1.0174 - This article can be downloaded from www.impactjournals.us 



96 



Jemilohun Vincent Gbenga 



27. K. Hinkelmann & O. Kempthorne, Design and analysis of experiments 1. (Second edition). 
(Wiley-Interscience (John Wiley & Sons), Hoboken, NJ. 2008). 

28. X. Wang, Weighted optimality of block design. Doctoral dissertation, Virginia Polytechnic Institute and 
State University, Virginia 2009. 



Index Copernicus Value: 3.0 - Articles can be sent to ecli(oi @ impact journals.us 



